ADP-01-44/T476 
FSU-CSIT-02-01 



Overlap Quark Propagator in Landau Gauge 



Frederic D.R. Bonnet 1 , Patrick 0. Bowman 2 , Derek B. Leinweber 1 , Anthony G. Williams 1 , 

J. B. Zhang 1 . 
1 CSSM Lattice Collaboration, 
Special Research Center for the Subatomic Structure of Matter ( CSSM) and Department of 

Physics and Mathematical Physics, University of Adelaide 5005, Australia 
2 Department of Physics and School for Computational Science and Information Technology, 
(N I Florida State University, Tallahasse FL 32306, USA 

(February 1, 2008) 

<\ 

o\ ' 

Abstract 

> : 

m . 

O ; 

The properties of the quark propagator in Landau gauge in quenched QCD 
are examined for the overlap quark action. The overlap quark action satisfies 
CN ■ the Ginsparg- Wilson relation and as such provides an exact lattice realiza- 

tion of chiral symmetry. This in turn implies that the quark action is free 
of 0(a) errors. We present results using the standard Wilson fermion kernel 
in the overlap formalism on a 12 3 x 24 lattice at a spacing of 0.125 fm. We 
| obtain the nonperturbative momentum-dependent wavefunction renormaliza- 

■ tion function Z(p) and the nonperturbative mass function M(p) for a variety 

of bare masses. We perform a simple extrapolation to the chiral limit for 
^ . these functions. We clearly observe the dynamically generated infrared mass 

and confirm the qualitative behavior found for the Landau gauge quark prop- 
agator in earlier studies. We attempt to extract the quark condensate from 
the asymptotic behavior of the mass function in the chiral limit. 



I. INTRODUCTION 

Hadron correlators on the lattice provide a direct means of calculating the physically 
observable properties of quantum chromodynamics (QCD). They are by construction color- 
singlet (i.e., gauge- invariant) quantities. Any finite, Boltzmann-distributed ensemble of 
gauge configurations has negligible probability of containing two gauge-equivalent configu- 
rations. Hence, there is a negligible probability of any gauge orbit being represented more 
than once in the Monte-Carlo estimate of the color-singlet hadron-correlator and this is the 
reason that there is no need to gauge fix in such calculations. 
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On the other hand, calculations of high-energy processes are carried out analytically 
with perturbative QCD, where it is necessary to select a gauge. Quark models and QCD- 
inspired Dyson-Schwinger equation models [|IJ are necessarily formulated in a particular 
gauge. The usual Fadeev-Popov gauge-fixing procedure is adequate for perturbative QCD. 
However in the nonperturbative infrared region standard gauge choices, such as Landau 
gauge, have Gribov copies, i.e., there are multiple gauge configurations on a given gauge 
orbit which satisfy the gauge-fixing condition. Since no finite ensemble will ever contain two 
configurations from the same gauge orbit, Landau gauge on the lattice actually corresponds 
to a gauge where there is a more or less random choice between the Landau gauge Gribov 
copies on the represented gauge orbits in the ensemble. Before gauge fixing, the ensemble 
contains configurations randomly located on their gauge orbits. After Landau gauge fixing 
on the lattice each configuration in the ensemble will typically settle on one of the nearby 
Landau gauge Gribov copies. This is the standard lattice implementation of Landau gauge 
and the one that we consider in this work. 

In order to study the transition from non-perturbative to the perturbative regime on the 
lattice we can study the gluon and quark propagators and vertices such as the quark- 
gluon vertex |1(J . By studying the momentum dependent quark mass function in the infrared 
region we can gain some insights into the mechanism of dynamical chiral symmetry breaking 
and the associated dynamical generation of mass. Studying the ultraviolet behaviour of 
propagators at large momentum is made difficult because of lattice artifacts causing the 
propagator to deviate strongly from its correct continuum behaviour in this regime. The 
method of tree-level correction was developed and used successfully in gluon propagator 
studies and has recently been extended to the case of the quark propagator !§-§[]. Some 
related studies have been performed for the case of domain- wall fermions J7|. Detailed 
discusions of nonpertubative renormalization for lattice operators can be found, e.g., in 
Refs. |i 

We present here results for the quark propagator obtained from the overlap quark action 
and using an improved gauge action and improved Landau gauge fixing. The overlap action 
is an exact realization of chiral symmetry on the lattice and is necessarily 0(a)-improved. 
In Sees. || and |TTT| we briefly introduce the improved gauge action and the lattice quark 
propagator respectively. In Sec. [IV| we introduce the overlap quark propagator and describe 
how it is calculated. Our numerical results are presented in Sec. [V] and finally in Sec. [VT| we 
give our summary and conclusions. 

II. IMPROVED GAUGE ACTION 

The tree-level O (a 2 ) -improved action is defined as, 

x [iv U x /i v 

where P^ u and R^ u are defined as 

P,u(x) = U,(x) U v {x + ft) UUx + 0)Ul(x), (2) 
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= U^(x) U v (x + p) U u {x + *> + //) Ufa + 1v) U ] u {x + v) Ul(x) 
+ UJx) U^x + fi) U v {x + 2ft) UUx + p. + u) UUx + u) U}(x). 



(3) 



The link product i? M „(x) denotes the rectangular 1x2 and 2x1 plaquettes. The mean link, 
Mo, is the tadpole (or mean- field) improvement factor that largely corrects for the quantum 
renormalization of the coefficient for the rectangles relative to the plaquette. We employ the 
plaquette measure for the mean link 



//,, = ( -Retr(P Ati/ (x)) 



1/4 



(4) 



where the angular brackets indicate averaging over x, «, and v. 

Gauge configurations are generated using the Cabbibo-Marinari [|n[] pseudo-heat-bath 
algorithm with three diagonal SU C (2) subgroups cycled twice. Simulations are performed 
using a parallel algorithm on a Sun Cluster composed of 40 nodes and on a Thinking Ma- 
chines Corporations (TMC) CM-5 both with appropriate link partitioning. We partition the 
link variables according to the algorithm described in Ref. [12]. We use 50 configurations 
generated on a 12 3 x 24 lattice at (3 = 4.60, selected after 5000 thermalization sweeps from a 
cold start and every 500 sweeps thereafter with a fixed mean-link value. Lattice parameters 
are summarized in Table |. The lattice spacing is determined from the static quark potential 
with a string tension yfa = 440 MeV [13|. 



The gauge field configurations are gauge fixed to the Landau gauge using a Conjugate 



Gradient Fourier Acceleration |L4| algorithm with an accuracy of 6 = \dfiA^(x)\ < 10 
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We use an improved gauge-fixing scheme to minimize gauge-fixing discretization errors. A 
discussion of the functional and method for improved Landau gauge fixing can be found in 
Ref. IBE 



III. QUARK PROPAGATOR ON THE LATTICE 

In a covariant gauge in the continuum the renormalized Euclidean space quark propagator 
must have the form 

{ ^ ,P) ii>A(C-p 2 ) + B(C,-p 2 ) ^^ + M(p 2 )' {) 

where ( is the renormalization point and where the renormalization point boundary condi- 
tions are Z((; ( 2 ) = 1 and M(( 2 ) = m(() and where m(() is the renormalized quark mass at 
the renormalization point. Since the gauge fixing condition has no preferred direction in color 
space, the quark propagator must be diagonal in color space, i.e., S^((; p) = S((; p)5^ where 
8 l i is the 3x3, SU(3) C identity matrix. The functions A((;p 2 ) and B((;p 2 ), or alternatively 

TABLE I. Lattice parameters. 



Action 


Volume 


A^Thcrm 


^Samp 


13 a (fm) 


n 


Physical Volume (fm 4 ) 


Improved 


12 3 x 24 


5000 


500 


4.60 0.125 


0.88888 


1.5 3 x 3.00 
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Z((;p 2 ) and M(p 2 ), contain all of the nonperturbative information of the quark propagator. 
Note that M(p 2 ) is renormalization point independent, i.e., since S(C;p) is multiplicatively 
renormalizable all of the renormalization-point dependence is carried by Z((;p 2 ). For suffi- 
ciently large momenta the effects of dynamical chiral symmetry breaking become negligible 
for nonzero current quark masses, i.e., for large ( and 7^ we have m(Q — > where 
is the usual current quark mass of perturbative QCD at the renormalization point £■ When 
all interactions for the quarks are turned off, i.e., when the gluon field vanishes, the quark 
propagator has its tree-level form 

5(0) G°) = ^^, (6) 

where m° is the bare quark mass. When the interactions with the gluon field are turned on 
we have 

5 (o)(p) _> S bare (a;p) = Z 2 (C;a)S(C;p) , (7) 

where a is the regularization parameter (i.e., the lattice spacing here) and Z 2 (C]a) is the 
quark wave-function renormalization constant chosen so as to ensure Z((;( 2 ) = 1. For 
simplicity of notation we suppress the a-dependence of the bare quantities. 

On the lattice we expect the bare quark propagators, in momentum space, to have a 
similar form as in the continuum except that the 0(4) invariance is replaced by a 

4-dimensional hypercubic symmetry on an isotropic lattice. Hence, the inverse lattice bare 
quark propagator takes the general form 



(8) 



We use periodic boundary conditions in the spatial directions and anti-periodic in the time 
direction. The discrete momentum values for a lattice of size Nf x N t , with rii = 1, .., iVj 
and n t = 1, .., N t , are 

2tt ( Ni\ , 2tt f 1 N t 



Defining the bare lattice quark propagator as 



(10) 



we perform a spinor and color trace to identify 



CM = ^tr^S^p)] and B(p) = ^tr^p)] . (11) 



v c 



The inverse propagator is 
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(S hare )- l (p) 



i(E M C M (p) 7l ,)+B(p) 



(12) 



C 2 (p)+S 2 (p) 

where C 2 (p) = E /t (^G°)) 2 - From Eq. fl|) we identify 



A. Tree-level correction 

At tree-level, i.e., when all the gauge links are set to the identity, the inverse bare lattice 
quark propagator becomes the tree-level version of Eq. (||) 

(S^)-\p) = i \^£C^(p) % j + B«0(p) . (14) 

We calculate (S^°'){p) directly by setting the links to unity in the coordinate space quark 
propagator and taking its Fourier transform 

It is then possible to identify the appropriate kinematic lattice momentum directly from 
the definition 

% = Cf{p) = (£(0)^)2% (B(o)(p))2 ■ ( 15 ) 

This is the starting point for the general approach to tree-level correction developed in earlier 
studies of the gluon propagator and the quark propagator ||-|(| . 

Having identified the appropriate kinematical lattice momentum q, we can now define 
the bare lattice propagator as 

^)= „.J +m = 1 j^ r WC-,*)S(C-, 1 ,) (16) 



and the lattice version of the renormalized propagator in Eq. 



' S ' ( ^> = ^ Lwr \ = 111 n • ( 17 ) 
,P) + B{(;p) i$ + M(p) 



The general approach to tree-level correction |2]|§-|6| utilizes the fact that QCD is asymp- 
totically free and so it is the difference of bare quantities from their tree-level form on the 
lattice that contains the best estimate of the nonperturbative information. For example, 
multiplicative tree-level correction for Z(p) and M(p) has the form 
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(18) 



The identification of the kinematical variable q ensures that A^°\p) = \/Z^°\p) = 1 by 
construction and so Z(jp) = Z^ c \p) and is already tree-level corrected. For overlap quarks 
we will see that M^°\p) = m° and so the mass function satisfies M(p) = M^ c \p) and needs 
no tree-level correction either. This feature is a major advantage of the overlap formalism. 



IV. OVERLAP FERMIONS. 



The overlap fermion formalism [[L6]-|25| realizes an exact chiral symmetry on the lattice 



and is automatically 0(a) improved, since any O(a) error would necessarily violate chiral 
symmetry | 20| . The massless coordinate-space overlap-Dirac operator can be written in 



dimensionless lattice units as [^T 



£(0) = ^[l+75#a] , (19) 

where H a is an Hermitian operator that depends on the background gauge field and has 
eigenvalues ±1. Any such D(0) is easily seen to satisfy the Ginsparg- Wilson relation [[27j 



{ 75 ,D(0)} = 2D(0) l5 D(0) (20) 

and, provided that its Fourier transform at low momenta is proportional to the momentum- 
space covariant derivative, it will satisfy a deformed lattice realization of chiral symmetry. 
It immediately follows from Eq. (|19D that 

Dt(0)D(0) = D(0)Dt(0) = i[Dt(0) + D(0)] (21) 

and that 

Dt(0)= 75^(0)75- (22) 

It also follows easily that {75, .D _1 (0)} = 275 and by defining .D _1 (0) = [.D _1 (0) — 1] we see 
that the Ginsparg- Wilson relation can also be expressed in the form 

{ l5 ,D-\0)} = 0. (23) 

The standard choice of H a (x,y) is H a = e(H w ) = H W /\H W \ = H w /(Hl J H w ) 1 ' 2 , where 
H w (x,y) = ^Dyj^x.y) is the Hermitian Wilson-Dirac operator and where D w is the usual 
Wilson-Dirac operator on the lattice. However, in the overlap formalism the Wilson mass 
parameter m w is the negative of what it is for standard Wilson fermions and at tree-level must 
satisfy < m w a < 2. In the overlap formalism m w is an intermediate lattice regularization 
parameter, it is not the bare quark mass. When interactions are present, we must have 
mi a < m w a in order that the Wilson operator has zero-crossings and, in turn, that D(0) 
has nontrivial topological charge. Numerical studies have found that mi ~ m c , where m, 



c 



is the usual critical mass for Wilson fermions PB1. The constraint m w a < 2 at tree level 
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arises from the fact that Wilson doublers reappear above this point. In summary, we use 
here H w (-m w ) = j 5 D w (-m w ). 

Recall that the standard Wilson-Dirac operator can be written as 

D w (x,y) = [(-m w a) + 4r] 8 x>y - ]- {(r - 7 M )^( X )^+A + ( r + 1t^) u l( x ~ ^)^y,x-fi} 



2k 



St. 



$x, y - K st 2j {(r - j^U^Sy^+p, + (r + 7 M )f/^(x - a/t)^-^} 



, (24) 



where the negative Wilson mass term (— m w a) is then defined by (-m K a) + 4r = l/2K st or 
equivalently 



2(-m w a) + (l//c c ) 



(25) 



and where k c throughout this work is the tree- level critical k, i.e., k c = l/(8r). 

In the present work we use the mean-field improved Wilson-Dirac operator, which can 
be written as 



D w (x,y) = 

2k 



UJx) 



C Z) 1 ( r - >)— W + (r + 7m) 
1 Wo 



£/t(x - a/2) 



"0 



J y,x~fi 



(26) 



We see that this is equivalent to the standard Wilson-Dirac operator with the identification 
of the mean-field improved coefficient k = K st UQ. It is U/uq that has a more convergent 
expansion around the identity than the links U themselves. The negative Wilson mass 
(—m w a) is then related to this improved k by 



K 



Wo 



2(-m w a) + (l//c c ) ' 



(27) 



The Wilson parameter is typically chosen to be r = 1 and we will also use r = 1 here in our 
numerical simulations. 

For this mean-field improved Wilson-Dirac choice we then have 



D(0) = -\l + D w (DlD w ) 



-1/2" 



(28) 



In coordinate space the Wilson-Dirac operator has the form D w = f+ (r/2)A + (—m w a), 
where V M is the symmetric dimensionless lattice finite difference operator, and A is the 
dimensionless lattice Laplacian operator. Recall that the Wilson mass term is (—m w a) here. 
Setting the links to the identity gives 



(l/2)(a A1 + d^hn + (r/2)(- d^d») + (-m w a) 



(29) 



where the partial derivatives are the forward and backward lattice finite difference operators. 
Hence we have from Eq. (|28[) that 
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f+ (r/2)A - m w a 



(30) 



y 



2m w a 



where the last line is a limit approached when operating on very smooth functions such that 
only first powers of derivatives are kept. The reason for needing a negative Wilson mass 
(—m w a) is now apparent, i.e., it is needed to cancel the 1 in -D(O). We see that at sufficiently 
fine lattice spacings and for pa <C 1 

D c (0) = (2m w )D(0) , (31) 

where D c (0) in the continuum limit becomes the usual fermion covariant derivative con- 
tracted with the 7-matrices, i.e., D c (0) — > Ip as a — > 0. 
The massless overlap quark propagator is given by 

s h -(o) ee d;\o) ee d;\o) - J- = J- [^(o) - 1] = ^-^(o) . (32) 

— II li uj —ill I'll) —ill I'fjj 

This definition of the massless overlap quark propagator follows from the overlap for- 
malism |T3 and ensures that the massless quark propagator anticommutes with 75, i.e., 
{75, S' bare (0)} = just as it does in the continuum [^l]]. At tree-level the momentum-space 



form of the massless propagator defines the kinematic lattice momentum q, i.e., we set the 
links to one such that we have for the momentum-space massless quark propagator 

S b ™ e (0,p) ee D-\0,p) - S(°\0,p) = I , (33) 

m 

where recall that p is the discrete lattice momentum defined in Eq. @ and q is the kinemat- 
ical lattice momentum defined in Eq. (p~5|) . We can obtain q numerically in this way from 
the tree-level massless quark propagator. We can compare this with the analytic form for q 
derived in the Appendix and given in Eq. ( |A9| ) . 

Note that for our mean-field improved Wilson-Dirac operator, the tree-level limit for 
defining q implies that we should take U — > / and uq — > 1 in D w while keeping k fixed. Thus 
the k that appears in the tree-level expression for g M in Eq. ( |A9|) is actually the improved 
k and not n st . This means that the tree- level Wilson mass parameter, m>w\ used in the 
Appendix is given by k = l/[2(— m^a) + (1/k c )] and hence differs from the m w in Eq. (p7|) 
used in the main body of the paper. We have found that the q obtained in this way gives a 
much superior large momentum behavior for the M(p) and Z(p) functions than is obtained 
when we do not use mean- field improvement. 

Having identified the massless quark propagator in Eq. (p2|), we can construct the massive 
overlap quark propagator by simply adding a bare mass to its inverse,i.e., 

Hence, at tree-level we have for the massive, momentum-space overlap quark propagator 
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S bare (m°,j9) -> S (0) (m°,p) 



(35) 



and the reason that the overlap quark propagator needs no tree-level correction beyond 
identifying q is now clear. 

In order to complete the discussion we now relate our presentation to standard notation 
used elsewhere. We first define the dimensionless overlap mass parameter 



m 
2m„ 



and then define D c 1 (fi) in analogy with Eq. ( |32| ) 

5 barc (m°) = D;\n) 



(36) 



(37) 



i.e., D~ 1 (n) is a generalization of D~ 1 (0) to the case of nonzero mass. Extending the analogy 
with the massless case we introduce D(n) and D(/j,), which are generalized versions of D(0) 
and -D(O), through the definitions 



D- 1 Qm) = ^D-\fi) and D^{p) = ^— [D~\n) - l] 



(38) 



We can now use Eq. ( [34]) and the above definitions to obtain an expression for D(fi). From 
Eqs. (151), (PO) and (EST) we see that we must have D(fi) — D(0) + fi. Inverting this gives 



1 



and so 

D-\p) = (l-fi) \£>(0)+H 
Inverting gives 

D(fi) = \D(0) + J \d(0) + 1 
and so finally 



[D-\ii)-1] = 6(0) +j* 



+ 1 



D(0) + 1 D(0)+fi 



-i 



and also D(0) = D(0) D(0) + 1 



n -1 



(39) 



(40) 



(41) 



D(fj)= (l-n)D(0)+n(D(0) + l) D(0) + 1 



[1 - fj)D{0) + /i=-[l + /i + (l- fj) l5 H a ] . 



(42) 



We have then recovered the standard expression for -D(/i) found for example in Ref. [21| and 
elsewhere. 

We see that the massless limit, m° — > 0, implies that /i — ► and D(fi) — > D(0), -D _1 (/i) — > 
.D _1 (0) and D~ l ([i) — > D ( T 1 (0). For non-negative bare mass m° we require // > 0. In order 
that the above expressions and manipulations be well-defined we must have \x < 1. Hence, 
< \i < 1 defines the allowable range of bare masses. 
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123456789 10 11 12 

p GeV 

FIG. 1. The kinematical momentum q for overlap quarks versus the discrete momentum p with 
both in GeV. No data cuts have been applied. The analytic result from the Appendix for the case 
of purely diagonal momenta is shown as the solid line for comparison. 

Our numerical calculation begins with an evaluation of the inverse of D(fi), where D(fi) 
is defined in Eq. ( |42"D and using H a = e(H w ) for each gauge configuration in the ensemble. 



We then calculate Eq. (|37|) for each configuration and take the ensemble average to obtain 
S hB,re (x, y). The discrete Fourier transform of this finally gives the momentum-space bare 
quark propagator, S ha,re (p), for the bare quark mass m°. 

Our calculations used k = 0.19 and u = 0.88888 and since at tree level k c = 1/8, we 
then have m w a = 1.661. Recall that the lattice spacing is a = 0.125 fm and so we have 
a -1 = 1.58 GeV and m w = 2.62 GeV. We calculated at ten quark masses specified by fi = 
0.024, 0.028, 0.032, 0.040, 0.048, 0.060, 0.080, 0.100, 0.120, and 0.140. This corresponds to 
bare masses in physical units of m° = 2fim w = 126, 147, 168, 210, 252, 315, 420, 524, 629, 
and 734 MeV respectively. 



V. NUMERICAL RESULTS 



We have numerically extracted the kinematical lattice momentum directly from the 
tree-level overlap propagator using Eqs. (|PT|), ( pT3| ) and (jl5|). In particular, by setting U — > / 



and Uq — > 1 we have numerically verified to high-precision the tree-level behavior in Eq. ( pB]) 
for all ten of our bare masses m°, which is a good test of our code for extracting the 

momentum-space quark propagator. We plot q = J ^2 q^ against the discrete lattice mo- 



mentum p = yJ^pPf,, in Fig. [1|. For pure Wilson quarks we would have g M = (1/a) sinfj^a). 

It is interesting that q for the overlap lies above the discrete lattice p, while q for Wilson 
quarks lies below. Of course in both cases, q — > p at small p. 
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123456789 10 11 12 

p GeV 

FIG. 2. The kinematical momentum q for overlap quarks versus the discrete momentum p with 
both in GeV. The cylinder cut has been applied and the hypercubic spread has been much reduced. 
The analytic result from the Appendix for the case of purely diagonal momenta is shown as the 
solid line for comparison. 

A. Data Cuts and Averaging 

To clean up the data and improve our ability to draw conclusions about continuum 
physics, we will on occasion employ the so-called "cylinder cut" , where we select only lattice 
four-momentum lying near the four-dimensional diagonal in order to minimize hypercubic 
lattice artifacts. This cut has been successfully used elsewhere in combination with tree-level 
correction in studies of the quark and gluon propagator |§f|||. It is motivated by the obser- 
vation that for a given momentum squared, (p 2 ), choosing the smallest momentum values 
of each of the Cartesian components, p M , should minimize finite lattice spacing artifacts. 

We calculate the distance Ap of a momentum four-vector p^ from the diagonal using 

Ap = \p\ sin 9 P , (43) 

where the angle 6 P is given by 

cos 6L = £ — -, (44) 

\v\ 

and n = |(1, 1, 1, 1) is the unit vector along the four-diagonal. For the cylinder cut employed 
in this study we neglect points more than one spatial momentum unit 2ir/Ni from the 
diagonal. To see that this cut has the desired effect of reducing hypercubic artifacts we 
plot the cut version of Fig. [I] in Fig. |2|. The cylinder-cut data points have much reduced 
hypercubic spread and lie on a smooth curve. We also sometimes make use of a "half-cut" 
where we only retain momentum components p^ half way out into the Brillouin zone. 

On an isotropic four dimensional lattice we have Z(4) invariance. Since our lattice is 
twice as long in the time-direction as it is in the spatial direction, this symmetry is broken 
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down to Z(3). This symmetry may be used to improve the statistics by averaging over Z(3)- 
identical momentum points. Since QCD and our lattice are parity invariant, we can also 
perform a reflection average at the same time. This average treats the negative momentum 
combinations in the same way as the positive ones. In an obvious notation, if we wish to 
calculate some quantity 5(1, 2, 3, 4), then we calculate all of the quantities 5(±i, ±j, ±k, ±4) 
for i, j, k being any permutation of 1, 2, 3 and perform an average over all of these quantities. 

We could also, in principle, average over all lattice starting points in the calculation of 
the propagator, since we should have translational invariance, i.e., S(y,x) should be the 
same for all equal (y — x). However, this averaging is too expensive to implement and we 
use a single starting point and calculate S(y, x) for all finishing points y^. We obtain the 
Fourier transform S(p) in the usual way with 

S(p)=J2^ ip<v ' x) S(y,x). (45) 



B. Overlap Quark Propagator 

In Fig. [3] we first show the half-cut results for all ten masses for both the mass and 
wavefunction renormalization functions, M(p) and Z^(p) = Z{C,\p) respectively, against 
the discrete lattice momentum p. Statistical uncertainties are estimated via a second-order, 
single-elimination jackknife. The renormalization point in Fig. |3] for Z^ K \p) has been chosen 
as C = 3.9 GeV on the p-scale. We see that both M(p) and Z^\p) are reasonably well- 
behaved up to 5 GeV although some anisotropy is evident. We see that at large momenta 
the quark masses are approaching their bare mass values as anticipated due to asymptotic 
freedom. 

In the plots of M{p) the data is ordered as one would expect by the values for /i, i.e., 
the larger the bare mass m° the higher is the M(p) curve. In the figure for Z^ K \p) the 
smaller the bare mass, the more pronounced is the dip at low momenta. Also, at small bare 
masses M(p) falls off more rapidly with momentum, which is understood from the fact that 
a larger proportion of the infrared mass is due to dynamical chiral symmetry breaking at 
small bare quark masses. This qualitative behavior is consistent with what is seen in Dyson- 
Schwinger based QCD models [|]] . The spread in the lattice data points indicates that some 
anisotropy from hypercubic lattice artifacts has survived the identification of the kinematical 
lattice momentum q. In Fig. [| we repeat these plots but now against the kinematical lattice 
momentum q. We see that the spread in the data is not significantly reduced and that the 
kinematical momentum reaches up to 12 GeV. 

The cylinder cut version of Figs. |3| and |] are given in Figs. |5] and ^| respectively. The 
cylinder cut removes almost all of the remaining spread in the overlap quark data and leaves 
data points which appear to lie on smooth curves. There is no apparent difference in the 
spread of the cylinder-cut data when plotted against p or q. Experience with the gluon 
propagator suggests that the continuum limit for Z(£,p) will be most rapidly approached 
as a — > by plotting it against its associated kinematical lattice momentum q. It is not 
obvious whether M(p) would converge to its continuum-limit behavior more rapidly as a — > 
by plotting against q or p or perhaps some other momentum scale. The only way to resolve 
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FIG. 3. The functions M(jp) and Z^\p) = Z(£;p) for renormalization point £ = 3.9 GeV (on 
the p-scale) for all ten bare quark masses for the half-cut data. Data are plotted versus the discrete 
momentum values defined in Eq. @), p = pfi, over the interval [0,5] GeV. The data in both 
parts of the figure correspond from bottom to top to increasing bare quark masses, i.e., \x = 0.024, 
0.028, 0.032, 0.040, 0.048, 0.060, 0.080, 0.100, 0.120, and 0.140, which in physical units corresponds 
to m° = 2fim w = 126, 147, 168, 210, 252, 315, 420, 524, 629, and 734 MeV respectively. The mass 
functions at large momenta are very similar to the bare quark masses as expected. 
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the g-scale) for all ten bare quark masses for the half-cut data. Data are plotted versus the discrete 
momentum values defined in Eq. (|T5|), q = q^, over the interval [0,12] GeV. The data in both 
parts of the figure correspond from bottom to top to increasing bare quark masses. The values of 
the bare quark masses are in the caption of Fig. |3[ 
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this question is to repeat the calculation on a lattice with different spacing a and to see 
which choice of momentum on the horizontal axis leads to a-independent behavior of Z(£;p) 
and M(p) most rapidly as a — > 0. This is left for future investigation. 



C. Extrapolation to Chiral Limit 

Our lightest bare quark mass is m° = 126 MeV and our heaviest is 734 MeV and hence 
we expect that our results should not be overly sensitive to the fact that our calculation 
is quenched. For exploratory purposes, we regard our simulation results at our bare quark 
masses to be a reasonable approximation to the infinite volume and continuum limits. We 
perform a simple linear extrapolation of our data to a zero bare quark mass, i.e., a linear 
extrapolation to the chiral limit. The results of our extrapolation for the mass function 
are shown in Fig. |7|. The top figure shows the chiral extrapolation of M(p) for the full 
uncut data set and plotted against p. The bottom figure shows the same results plotted 
against q up to 12 GeV. The fact that the linear extrapolation gives an M(p) which vanishes 
within statistical errors at large momenta confirms that our simple linear extrapolation is 
reasonable at large momenta. In fact, the data are found to be consistent with such a linear 
fit at all momenta for the bare masses considered. 

In Fig. |8] we plot the cylinder-cut data after the linear chiral extrapolation for both 
functions M(p) and Z^ R \p) = Z{C\p). These are shown against both p and q with the 
renormalization points chosen as in the previous figures, i.e., 3.9 GeV and 8.2 GeV for plots 
against p and q respectively. We see that both M(p) and Z^(p) deviate strongly from the 
tree-level behavior. In particular, as in earlier studies of the Landau gauge quark propagator 
0-0, we find a clear signal of dynamical mass generation and a significant infrared suppres- 
sion of the Z((;p) function. At the most infrared point, the dynamically generated mass has 
the value Mir, = 297(11) MeV and the momentum-dependent wavefunction renormalization 
function has the value Zir = 0.48(2). The result Mir = 297(11) MeV is similar to typical 
mass values attributed to the "constituent quark" mass and is approximately 1/3 of the 
proton mass. These values are very similar to the results found in previous studies 
and are also similar to typical values in QCD-inspired Dyson-Schwinger equation models 
0JHJ1 



As the bare mass m° is increased from the chiral limit as in Figs. [5] and |6|, we are increasing 
the proportion of explicit to dynamical chiral symmetry breaking. Associated with this we 
see that the dip in Z((; p) becomes less pronounced and the relative importance of dynamical 
chiral symmetry breaking in the mass function M(p) also decreases, i.e., we see that M(p) 
becomes an increasingly flat function of momentum as the bare mass is increased. 

In the continuum at one loop order in perturbation theory and in the presence of ex- 
plicit chiral symmetry breaking (i.e., a nonzero bare mass), the asymptotic behavior of the 
mass function is that of the running quark mass. Specifically for large, Euclidean p 2 and 
renormalization point ( we have the one-loop result [0 



M(p 2 ) P M°° m c 



ln(C 2 /A 2 



QCDy 



ln(p 2 /A^c D ) 



(46) 
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FIG. 5. The functions M{jp) and Z^\p) = Z(£;p) for renormalization point £ = 3.9 GeV (on 
the p-scale) for all ten bare quark masses and for data with a cylinder-cut, i.e., the data is identical 
to that of Fig. || except that it has been cylinder cut (one spatial momentum unit) rather than 
half-cut. 
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FIG. 7. The chiral limit mass function M(p) obtained from a simple linearly extrapolation of 
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where dyi = 12/(33 — 2Nf) is the anomalous mass dimension, is the current quark mass, 
Nf is the number of quark flavours, and Aqcd is the QCD scale parameter.. In this limit we 
see then that the mass at the renormalization point, m(() = M(( 2 ), approaches the current 
quark mass, i.e., m(() = M(( 2 ) — » as stated earlier. The vanishing of the bare mass 
m° defines the chiral limit and in that case the current quark mass also vanishes and the 
asymptotic behavior of the mass function at one-loop becomes 

3 [ H ( 2 /A 2 QCD )] d "P* 1 Q J 

This is the asymptotic behavior of the dynamically generated quark mass. We see that the 
running mass in Eq. (^6|) falls off logarithmically with momentum, whereas from Eq. ( f4"?|) 
the dynamically generated mass falls off more rapidly (as 1/p 2 up to logarithms) in the 
chiral limit. This is the reason that the effects of dynamical chiral symmetry breaking can 
be neglected at large momenta. Since M(p 2 ) is renormalization point independent, the 
combinations m^[ln(^ 2 /AQ CD )] dM , (qq)(/ [ln(£ 2 /AQ CD )] M , and m^(qq)^ are renormalization 
group invariant. The anomalous dimension of the quark propagator itself vanishes in Landau 
gauge. Hence in the continuum in Landau gauge 

Z((;p 2 ) p2 ^° 1. (48) 

In our lattice results we clearly observe that Z^\p) = Z[C,\ p) behaves in a way consistent 
with Eq. (|48"D . We can then also examine whether or not the asymptotic behavior of our 
linearly extrapolated chiral mass function satisfies Eq. ([|7]). Since we are working in the 
quenched approximation we have Nf = 0. We attempt to extrapolate the quark condensate 
for three different values of Aqcd, i.e., 200, 234, and 300 MeV. Note that 234 MeV is among 
typical values quoted for quenched QCD p3| . 

We also do the extraction over different ultraviolet fitting windows in order to verify the 
insensitivity of the chiral condensate to the fitting window. Since it is at present unclear 
whether M(p) most rapidly approaches the continuum limit when plotted against p or plotted 
against q, we have performed the fit to both, i.e., we have fitted Eq. ([4^) to both sets of 
ultraviolet data for M(p) in the half-cut version of the data in Fig. |7]. 

A summary of the fitting results is shown in Tables |T| and [TTJ for various fitting regions 
and QCD scale parameters. As is standard practice, we quote the extracted condensate 
at the renormalization scale ( = 1 GeV using the renormalization scale independence of 
{qq)(/ \\n(( 2 / Aq CD )Y m ■ It is the latter that is extracted in the fit to the chiral M{p). The 
extracted condensate is relatively insensitive to the value of Aqcd and the fitting window. 
There is however a very strong dependence on which momentum scale is used for M(p), i.e., 
~ 350 MeV for p compared with ~ 600 MeV for q. It is clear that a quantitatively meaningful 
extraction of the quark condensate will require us to establish which momentum scale for 
M{p) most rapidly reproduces the continuum limit as a — ► 0. Other attempts [PT[-|3"3"| to 
directly calculate the quark condensate in the overlap formalism in quenched QCD suggest 
a condensate value ~ 250 MeV, which implies that M(p) may be more appropriately plotted 
against the discrete lattice momentum p. This resolution of this issue is beyond the scope of 
the present study and is left for future work. However, once the correct momentum scale is 
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TABLE II. Summary of the results for the quark condensate, — (qq)^ , extracted from Eq. ( |47D 
in MeV and scaled to the renormalization point £ = 1.0 GeV. The fit was done using Eq. ( |47| ) for 
each of Aqcd = 200, 234, 300 MeV on various momentum windows using M(p) plotted against the 
discrete lattice momentum p. 



p GeV 


200 MeV 


Aqcd 
234 MeV 


300 MeV 


3- 5 

4- 5 


356(35) 
352(69) 


347(34) 
344(67) 


333(32) 
330(64) 


1/3 

TABLE III. Summary of the results for the quark condensate, — (qq)^ , extracted from 
Eq. (f4~7|) in MeV and scaled to the renormalization point £ = 1.0 GeV. The fit was done us- 
ing Eq. (|47|) for each of Aqcd = 200, 234, 300 MeV on various momentum windows using the 
linearly extrapolated half-cut data for M(p) plotted against the kinematical lattice momentum q. 


p GeV 


200 MeV 


Aqcd 
234 MeV 


300 MeV 


3-5 
3-7 
3-9 


604(68) 
600(66) 
594(63) 


591(67) 
587(65) 
581(62) 


566(64) 
562(62) 
557(59) 


4-5 
4-7 
4-9 


613(81) 
600(71) 
589(67) 


599(79) 
586(70) 
576(65) 


574(76) 
563(67) 
553(63) 


5-7 
5-9 


590(66) 
577(41) 


577(65) 
563(62) 


556(32) 
541(60) 



identified and the continuum limit estimated, the good quality of the overlap data indicate 
that an extraction of the quark condensate should be possible. 

VI. SUMMARY AND CONCLUSIONS 

To the best of our knowledge, this is the first detailed study of the Landau gauge 
momentum-space quark propagator in the overlap formalism. By construction, the over- 
lap quark propagator needs no tree-level correction beyond the identification of the appro- 
priate kinematical lattice momentum q. The quality of the data in the overlap formalism 
is seen to be far superior to that from earlier studies which use an C(a)-improved 

Sheikholeslami-Wohlert (SW) quark action with a tree-level mean-field improved clover co- 
efficient c sw . In these earlier studies it was found that careful tree-level correction schemes 
are essential and the resulting corrected data remain of inferior quality. The quality of the 
data for the improved staggered quark action, the so-called 'Asqtad" action with 0(a 4 ) and 
0(a 2 g 4 ) errors, is also seen to be superior to the C(a)-improved quark action and these 
results are qualitatively consistent with what we have found here. 



21 



We use ten different quark masses and observe an approximately linear relation between 
the bare mass and the current quark mass for bare masses in the range ~ 125 to ~ 730 MeV. 
This allows a simple linear extrapolation to the chiral limit. Such a linear extrapolation is 
justified in the ultraviolet, since the resulting ultraviolet mass function vanishes within errors 
in the chiral limit as expected. For the most infrared momentum points in the chiral limit 
using this linear extrapolation we find Mir = 297(11) MeV and = 0.48(2) for the mass 
function and the momentum-dependent wavefunction renormalization function respectively. 

An extraction of the quark condensate from the ultraviolet behavior of the chiral extrap- 
olated mass function is possible with this quality of data. However, it is clear that this can 
not be done in a quantitatively reliable way until one or more additional lattice spacings 
become available so that we can identify the appropriate momentum against which to plot 
M(p 2 ). 

The first calculation presented here is performed on a relatively small volume of 1.5 3 x 
3.0 fm 4 and an intermediate lattice spacing of 0.125 fm. Ultimately, a variety of lattice 
spacings and volumes should be used so that a study of the infinite- volume, continuum limit 
of the overlap quark propagator can be performed. It will also be interesting to simulate at 
lighter quark masses in order to study the chiral limit of the quenched theory in some detail. 
Finally, one should consider kernels in the overlap formalism other than the pure Wilson 



kernel, e.g., using a fat-link irrelevant clover (FLIC) action [34] as the overlap kernel [35 
These studies are currently underway and will be reported elsewhere. 
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APPENDIX A: TREE-LEVEL BEHAVIOR 

1. Tree-level overlap propagator 

We can derive an explicit form for the tree- level (i.e., free) overlap quark propagator 
with the Wilson fermion kernel. Let us work in the infinite volume limit with finite lattice 
spacing a. 

It is convenient to define the dimensionless momentum variables 

k^ = sm(p^a) , kf, = 2 sin(p M a/2) . (Al) 

Let us also define the dimensionless combination 

(A2) 



A 



-am 



(0)' 



2 . 
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where (—am$) is the negative, dimensionless tree- level Wilson mass defined by n = 
l/[2(— rriw^a) + (1/k c )]. We have k c = 1/8 and r is the usual Wilson parameter, (typi- 
cally one chooses r = 1). Note that for small momenta we have A < . We can then write 
the momentum-space Wilson operator at tree-level as 



1 

2k 



D„. = — r -k ■ .l) 



It follows that 



dLd,. 



2k 



(A3) 



(A4) 



where it is to be understood that by definition only the positive root is kept. In Euclidean 
space it is clear that we will always have k + A 2 > and the square root is always well- 
defined. The momentum-space overlap Dirac operator can then be written as 



D(0) 



1 + 



D„ 



VdId v 



1 + 



ry • k + A 



Vk 2 + A 2 ^ 
i7- k+ {A+ V~k 2 + A 2 } 

V~k 2 + A 2 



(A5) 



We see that H a = ^D w j a/ D\,D w and that H\ = H a and H 2 = 1 as required in the overlap 
formalism, i.e., H a has eigenvalues ± 1 . We can readily invert D(0) to give 



D-\0) = 2Vk 2 + A 2 
— ry • k 



A+V~k 2 + A 2 



-i7- jfe+ {A+ ^k 2 + A 2 } 
k 2 + {A+y/k 2 + A 2 Y 

+ 1 



(A6) 



and then 



D _1 (0) = D-\0) - 1 = 



-ry ■ k 



A + V~k 2 + A 2 « 7 • k |A + Vk 2 + A 2 } ' 



(A7) 



Clearly then {D 1 (0),75} = as it must in the overlap formalism. The tree- level 
momentum-space overlap quark propagator in the massless limit is then given by 



2ml 0) 



2mi 0) 



(A8) 
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Hence, we recognize that the kinematical tree-level momentum is given by 

\A+ Vk 2 + A 2 } 



q, = 2mfk^ j- >- . (A9) 

k 



This analytic form for the q = y^2q 2 versus p = yY^Pft behavior is plotted as the solid 

line in Figs, [l] and |2| for the case of purely diagonal momenta, (p± = P2 = P3 = Pa)- The 
analytic form can be checked against each (diagonal or non-diagonal) point on a case-by-case 
basis and they agree to within numerical precision. 

We can verify analytically that g M — > p^ as p — > as seen in Fig. |2] . Note that as p — ► 
we have A < and hence 



k 2 + \A\ 2 -\A\ 
P 



(0) - 

^7Ti y _ sin(p M a) p^ . (A10) 



2. Tree-level dispersion relation 



The massless, tree-level overlap propagator has the momentum-space form 



S (0) (0,p) 



(All) 



and so has poles when q 2 = . We can analytically continue q^ — > iE and then we find poles 
at E = \q\, i.e., in terms of our tree- level corrected propagator we have a perfect massless 
dispersion relation. 

However, for hadronic properties without tree-level correction it is the dispersion relation 
in p that is relevant, i.e., we need to analytically continue p^ — > iE and find the poles in 
S^(0,p). Our discussion here generalizes that given in Ref. [20|. It is clear from Eq. (|A4|) 
and Eq. ( |A~5| ) that the analytic continuation is only defined in the region where k 2 + A 2 > 0, 
since otherwise the argument of the square- root is negative and the definition of -D(O) has 
no meaning. 

The poles occur when 



= q 2 = 4(m 



(Oh 2 



|A+ \fk 2 



+ A 2 



4rf) 2 (l + =f 



2A 2 



k 2 



l + sgn(A)Jl + (k 2 /A 2 



(A12) 



Provided A < we see that q 2 — > as k 2 — > 0. Consider these poles when p 4 — > iE and 
p = (0, 0,p), then the conditions k = , A < become 
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sin 2 (iEa) + sin 2 (pa) = , (A13) 

(o) 



sin 2 (z£a/2) + sin 2 (pa/2) < (A14) 
respectively. Thus we have poles at 

cosh(Ea) = y^l + sin 2 (pa) (A15) 

when we satisfy the condition 

(o) 

cosh(Ea) > 2 - cos(pa) . (A16) 

r 

Note that the analytic continuation to Minkowski space is only well-defined when the square- 
root operation is well-defined, i.e., for k 2 +A 2 > and so this condition must also be satisfied. 
We can rewrite this condition as 

2 [2 — (ami} ) — cos(pa)] 



Recall that Eq. ( |A16|) is equivalent to the condition A < 0. If A = then q 2 = 4(m£^) 2 ^ 
for any real k 2 and hence there are no poles. If A > then in the region where the square-root 

is well defined |yl + + A 2 1 > and there are no poles in that case either. 
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FIG. 9. The dispersion relation for the overlap quark propagator of Eq. ( |A15| ) is shown as 
the solid line and corresponds to k = 0. The dispersion relation does not continue into the region 
where A > 0, i.e., it does not extend to the right beyond the long-dash dot line denoting ^4 = 0, 



[the solution of Eq. (A16)]. The analytic continuation to Minkowski space has no meaning when 
k 2 + A 2 < 0, i.e., it is undefined above the short-dash dot line, [i.e., the solution of Eq. (A17)]. 
The intersection point of these three curves is where we simultaneously have A, k , and k + A 2 
equal to zero. Also shown for reference are the dispersion relations for the continuum limit (short 
dashes) and for the ordinary Wilson action (long dashes). 
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